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O 

^\ ■ We prove an existence and uniqueness result for the obstacle 

problem of quasilinear parabolic stochastic PDEs. The method is 
s ! . based on the probabilistic interpretation of the solution by using the 

backward doubly stochastic differential equation. 

^ ■ 1. Introduction. We consider the following stochastic PDE, in 

^ ■ dut(x) + {kAu t (x) + f t (x,ut(x), Vu t (x)) 

(1) +div g t (x,Ut(x),Vu t (x))]dt 

+ h t (x,Ut(x),Vu t (x)) ■ tl3 t = 0, 



> 



id 



over the time interval [0, T], with a given final condition ut = $ and f,g = 
(51 j • • • 1 9d)> h = (h\, . . . , h d i) nonlinear random functions. The differential 
term with dBt refers to the backward stochastic integral with respect to a 
d 1 -dimensional Brownian motion on (Q, T , P, (Bt)t>o)- We use the backward 
notation because in the proof we will employ the doubly stochastic frame- 
work introduced by Pardoux and Peng [16] (see also Bally and Matoussi [2] 
and Matoussi and Xu [13]). 

In the case where / and g do not depend of u and Vu, and if h is identically 
null, the equation (1) becomes a linear parabolic equation, 

'. (2) d t u(t,x) + ±Au(t,x) + f(t,x)+divg(t,x) = 0. 
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If v. [0,T] x M — > M. is a given function such that v(T,x) < we may 

roughly say that the solution of the obstacle problem for (2) is a function u 6 
L 2 ([0, T]; such that the following conditions are satisfied in (0,T) x 

R d : 



The relation (ii) means that the distribution appearing in the LHS of the 
inequality is a nonpositive measure. The relation (iii) is not rigourously 
stated. We may roughly say that one has dtu + \ Au + / + div g = on the 
set {u > v}. 

If one expresses the obstacle problem for (2) in terms of variational in- 
equalities one should also ask that the solution has a minimality property 
(see Bensoussan-Lions [3], page 250, or Mignot-Puel [14]). 

The work of El Karoui et al. [9] treats the obstacle problem for (2) within 
the framework of backward stochastic differential equations (BSDE in short). 
Namely, the equation (2) is considered with / depending of u and Vu, while 
the function g is null (as well h) and the obstacle v is continuous. The so- 
lution is represented stochastically as a process and the main new object 
of this BSDE framework is a continuous increasing process that controls 
the set {u = v}. This increasing process determines in fact the measure 
from the relation (ii). Bally et al. [1] point out that the continuity of this 
process allows one to extend the classical notion of strong variational so- 
lution (see Theorem 2.2 of [3], page 238) and express the solution to the 
obstacle as a pair (u, v) where v equals the LHS of (ii) and is supported 
by the set {u = v }. Moreover, based on this observation Matoussi and Xu 
[12] generalized the work under monotonicity and general growth conditions. 
They have also used the penalization method and stochastic flow technics 
(see [2] and [11] for more details on this method). In the present paper, 
we similarly consider the solution as a pair (u,u), point of view which 
has the advantage of expressing the notion of solution independently of 
the double stochastic framework and without the minimality property of 
Mignot-Puel [14], which would be very difficult to manipulate in the case 
of the stochastic PDE. In Section 2.2, we are going to examine the po- 
tential and the measure associated to a continuous increasing process. We 
call such potentials and measures, regular potentials, respectively regular 
measures. 

Now let us consider the final condition to be a fixed function $ £ L 2 (R rf ) 
and the obstacle v be a random continuous function, u:!!x [0, T] x M. d — > R. 



(3) 



(i) 
(ii) 

(iii) 
(iv) 



u>v, dt(E) dx-a.e., 
d t u + i An + / + div g < 0, 

(u — v)(dtu + ^An + / + div g) 
ut = dx-a.e. 
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(4) 



Then the obstacle problem for the equation (1) is defined as a pair (u, u), 
where v is a random regular measure and u G L 2 (f2 x [0, T];H l (M. d )) satisfies 
the following relations: 

(i') u > v, dP ® dt <g> dx-a.e., 

(ii') du t (x) + [^Ant(x) + f t (x,u t (x),Vu t (x)) 

+ divg t (x,u t (x),Vu t (x))] dt 

+ h t (x, Ut(x),Vut(x)) ■ cLB t = —v(dt, dx) a.s., 

(in') v(u > v) = a.s., 

(iv') ut = ®, dF® dx-a.e. 

In Section 2.4, we explain the rigorous sense of the relation (iii') which is 
based on the quasi-continuity of u. The main result of our paper is Theorem 
4 which ensures the existence and uniqueness of the solution of the obstacle 
problem for (1). The method of proof is based on the penalization procedure 
and the doubly stochastic calculus which is essential, although the definition 
of the solution and the statement of the result avoids the doubly stochastic 
framework. 

Similarly to the case treated in El Karoui et al. [9], the most difficult 
point is to show that the approximating sequence converges uniformly on the 
trajectories over the coincidence set {u = v}. This is proven in Lemma 7. The 
existence and uniqueness of the solution for equation (1) (without obstacle) 
has already been proven in [7]. An essential ingredient in the treatment of 
the quasilinear part is the probabilistic representation of the divergence term 
obtained in [17] as well as the doubly stochastic representation corresponding 
to the divergence term of the stochastic PDE in [7]. We must mention the 
work of Nualard and Pardoux [15] and Donati-Martin and Pardoux [8] who 
studied a particular class of obstacle problem for stochastic PDE driven by 
some space-time white noise by using a different techniques. 

Finally, we would like to thank our friend Vlad Bally for a stimulating 
discussion on the obstacle problem we had "la Gare de Montparnasse" and 
the referee for helping us to improve the presentation. 

2. Preliminaries. The basic Hilbert space of our framework is L 2 (M' i ), 
and we employ the usual notation for its scalar product and its norm, 

. 1/2 

(u,v) = I u(x)v(x)dx, IMh — ( / u 2 (x)dx 
In general, we shall use the notation 

(u,v)= / u(x)v(x) dx, 
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where u, v are measurable functions denned in M. d and uv £ L 1 (M d ). 

Our evolution problem will be considered over a fixed time interval [0, T] 
and the norm for a function L 2 ([0,T] x M. d ) will be denoted by 

IM|2,2=( / / \u(t,x)\ 2 dxdt) 
\Jo JR d J 

Another Hilbert space that we use is the first order Sobolev space H l (W i ) = 
H^MJ^). Its natural scalar product and norm are 

(u,v) H i {Rd) = (u,v) + (Vu,Vv), IMItfi(Rd) = (H|! + l|Vn||^) 1/2 , 

where we denote the gradient by Vu(t,x) = (d\u(t,x), . . . ,ddu(t,x)). 

Of special interest is the subspace F C L 2 ([0,T];// 1 (]R c( )) consisting of 
all functions u(t,x) such that t^-ut = u(t,-) is continuous in L 2 (M rf ). The 
natural norm on F is 

IM|t= sup ||n t ||2+( / \\Vu t \\ 2 dt 

0<t<T \Jo 

The Lebesgue measure in R d will be sometimes denoted by m. The space 
of test functions which we employ in the definition of weak solutions of the 
evolution equations (1) or (2) is V T = C°°([0, T]) <g> C™(R d ), where C°°([0,T]) 
denotes the space of real functions which can be extended as infinite differ- 
entiable functions in the neighborhood of [0,T] and is the space of 

infinite differentiable functions with compact support in M. d . 




2.1. The probabilistic interpretation of the divergence term. The opera- 
tor dt + ^A, which represents the main linear part in the equation (1), is 
probabilistically interpreted by the Brownian motion in M. d . We shall view 
the Brownian motion as a Markov process, and therefore we next introduce 
some detailed notation for it. The sample space is W = C([0, oo); Mr), the 
canonical process (Wt)t>o is defined by Wt(ui) = u(t), for any u £ SI', t > 
and the shift operator, 9 t : Q' — > fl\ is defined by 6t(uj)(s) = co(t + s), for any 
s > and t > 0. The canonical filtration = a(W s ;s < t) is completed by 
the standard procedure with respect to the probability measures produced 
by the transition function 

P t (x, dy) = q t (x -y)dy, t > 0, x G R d , 

where qt(x) = {2'Kt)~ d / 2 exp(— \x\ 2 /2t) is the Gaussian density. Thus, we get 
a continuous Hunt process (fi', Wt,6 t ,F ,Ft,^ x ). We shall also use the back- 
ward filtration of the future events F[ = a(W s ;s > t) for t > 0. P° is the 
Wiener measure, which is supported by the set Q' = {u E £l',w(0) = 0}. 
We also set IIo(w)(t) = co(t) — co(0),t > 0, which defines a map IIo : — > fl' . 
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Then II = (W , LT ) : Q' R d x % is a bijection. For each probability mea- 
sure on M d , the probability P M of the Brownian motion started with the 
initial distribution \i is given by 

p^ = n- 1 (/i^p°). 

In particular, for the Lebesgue measure in M. d , which we denote by m = dx, 
we have 

F m = Il~ 1 (dx®F ). 

These relations are saying that Wq is independent of Ilo. It is known that 
each component (Wf)t>o of the Brownian motion, i = 1, . . . , d, is a martin- 
gale under any of the measures P^. The next lemma shows that (Wf_ r ,J-l_ r ), 
r 6 (0,t], is a backward local martingale under P m . 

Lemma 1. Let < s < t. If A £ a(W t ) is such that E m [\W t \;A] < oo, 
then one has E m [|W s |; A] < oo. Moreover, for each B £ T[, and i = 1, . . . ,d, 
one has 

E m [W l s ; AnB]= E m [Wp, Af]B]. 

Proof. We note that W% is uniformly distributed, and consequently for 
each c > 0, the set A c = {\Wt\ < c} satisfies 

E m [|H^M c ]<oo. 

This shows that the class of the sets to which applies the statement is rather 
large. 

The vector (Wo, W s - W Q ,W t - W s ) has the distribution m <g) jV(0,s) ® 
M(0,t — s), under the measure P m . Then one deduce that (W s ,W t — W s ) 
has the distribution m <S> Af(0,t — s) and we may write, for <fi,(f2 £ C c (R d ), 

E m [i Pl (Wt-W s )ip 2 (W t )}= [ [ My)Mx + y)qt-s(y)dydx 

= (/ ^2{x)dx)[ I yi(y)q t - s (y)dy 
\JR d J \J R d 

This relation shows that the vector (Wt — W s ,Wt) has the distribution 
W(0,i — s) <8>m, under P m .Then the obvious inequality |W S | < |Wt| + |Wf — 
W a |(l/|w t i<i} + \Wt\) allows one to deduce the first assertion of the lemma. 
In order to check the second assertion of the lemma, we write 

E m [W*;A n B] = E m [Wi; AnB}- E m [W l t - W l s ; Af]B] 

and all that it remains to check is that the last term is null. In order to show 
this, one first observes that the distribution of the vector (Wt — W s ,Wt, W t i — 
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W t , W t 2 - W t i,...,W t n - W t n-i) is M(0,t - s) m A/"^ 1 - t) ••• 
7V(0, i n — t n_1 ), for each system s < t < t 1 < • • • < t n . Then one has, for each 
B e a(W t i - W t , . . . , W t n - W tn -i), 

E m [Wt - Wl; A(~)B]= E°[Wi - Wi]m(A)T° (B) = 0, 

which implies the assertion of the lemma. □ 

Now let us assume that / and \g\ belong to L 2 ([0,T] x M. d ) and u € F is 
a solution of the deterministic equation (2). Let us denote by 

(5) f g r *dW r = J2( f 9i(r, W r ) dW l r + f 9i (r, W r ) 0W) 

Then one has the following representation (Theorem 3.2 in [17]). 

Theorem 1. The following relation holds F m -a.s. for each < s < t < 

T: 

(6) u t {Wt)-u s (W s )=Y J d lUr (W r )dW l r - f r {W r )dr-- g r *dW r . 

■ i i/ S J S 'J s 



1=1 



In [17], one uses the backward martingale < M^>' 1 defined under an arbitrary 
P^, with [i a probability measure in M. d , in order to express the integral 
f t g r * dW r . Though formally the definition looks different, one easily sees 
that it is the same object. 

2.2. Regular measures. In this section, we shall be concerned with some 
facts related to the time-space Brownian motion, with the state space [0, T[xM a! , 
corresponding to the generator dt + |A. Its associated semigroup will be de- 
noted by (Pt)t>o- We may express it in terms of the Gaussian density of the 
semigroup (Pt)t>o m the following way: 



P t tp(s,x) 



/ Qt(x,y)i>(s + t,y)dy, if s + t<T, 
, 0, otherwise, 



where tp : [0, T fxPJ^Misa bounded Borel measurable function, s £ [0, T[, x € 
M. d and t > 0. So we may also write (Pttp) s = Ptipt+s if s + t <T. The cor- 
responding resolvent has a density expressed in terms of the density qt too, 
as follows: 

U a if)(t,x)= [ [ e- a(s - t) q s -t(x-y)4>(s,y)dyds 
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or 



T 



In particular, this ensures that the excessive functions with respect to the 
time-space Brownian motion are lower semicontinuous. In fact, we will not 
use directly the time space process, but only its semigroup and resolvent. 
For related facts concerning excessive functions, the reader is referred to [4] 
or [6]. Some further properties of this semigroup are presented in the next 
lemma. 

Lemma 2. The semigroup (Pt)t>o acts as a strongly continuous semi- 
group of contractions on the spaces L 2 ([0, T[xR d ) = L 2 ([0, T[; L 2 (R d )) and 
V(iO,T[;H\R d )). 

Proof. Obviously, it is enough to check the following relations: 



First, we note that for each function u E L 2 ([0, T[xM rf ) and r > 0, one has 



This property is obvious for a function u £ C c ([0,T[xR ) and then it is 
obtained by approximation for any function in L 2 ([0, T[xM rf ). Then the re- 
lation 



easily follows. From it, one deduces the strong continuity of (Pt)t>o on 
L 2 ([0,r[xR d ). 

In order to prove the same property in the space L 2 ([0,T[;ff 1 (M c( )), one 
should start with the relation 



which holds for each u £ C£°([0, T[xM rf ) and then repeat, with obvious mod- 
ifications, the previous reasoning. □ 

The next definition restricts our attention to potentials belonging to F, 
which is the class of potentials appearing in our parabolic case of the obstacle 
problem. 
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Definition 1. (i) A function ip : [0,T] x R d — > R is called quasicontinu- 
ous provided that for each e > 0, there exists an open set, D £ C [0, T] x R rf , 
such that if) is finite and continuous on D c e and 

F m ({uj£n'\3t€ [0,T] s.t. (t,Wt(w)) £D £ }) < e. 

(ii) A function u : [0, T] x R d — > [0, oo] is called a regular potential, provided 
that its restriction to [0,T[xR d is excessive with respect to the time-space 
semigroup, it is quasicontinuous, u E F and \mit^,T u t = in L 2 ' 



Observe that if a function ip is quasicontinuous, then the process 
{'4 , t{W't))te[o,T\ is continuous. Next, we will present the basic properties of the 
regular potentials. Do to the expression of the semigroup (Pt)t>o i n terms of 
the density, it follows that two excessive functions which represent the same 
element in F should coincide. 

Theorem 2. Let u£ F. Then u has a version which is a regular poten- 
tial if and only if there exists a continuous increasing process A = (^)te[o,ri 
which is (Ft)t£\o,T]' a d a ptsd and such that Aq = 0, E m [A|,] < oo and 

(i) u t {W t ) = E[A T \F t ] - A t F m -a.s. 

for each t G [0, T]. The process A is uniquely determined by these properties. 
Moreover, the following relations hold: 



d ,. T 

(ii) u t (W t ) =A T -A t -J2 d iUs {W s ) dWt 

i=i 

(hi) Hutlll + £ \\Vu s \\l ds = E m (A T - A t ) 2 , 
(iv) (u ,(p ) + J \^(Vu s ,V<p s ) + {u s ,d s (p s )\ds 



im -a.s., 



T 







ip(s, x)u{ds dx) 



for each test function (p S T> T , where v is the measure defined by 

(v) v(ip)=E m [ T p(t,Wt)dA t , ^C c ([0,T]xR d ). 



o 



Proof. We first remark that the uniqueness of the increasing process 
in the representation (i) follows from the uniqueness in the Doob-Meyer 
decomposition. 
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Let us now assume that u is a regular potential which is a version of u. 
We will use an approximation of u constructed with the resolvent. By the 
resolvent equation, one has 

all a u = alio (u — aU a u) . 

Let us set f n = n(u — nU n ll) and u n = nU n u = Uof n . Since u is excessive, 
one has f n > and u n ,n £ N*, is an increasing sequence of excessive func- 
tions with limit u. In fact u n ,n £ N*, are potentials and their trajectories 
are continuous. On the other hand, the trajectories t — > ut(Wt) are contin- 
uous on [0,T[ by the quasi-continuity of u. The process {ut(Wt))t€.[Q,T[ is a 
super-martingale, and because \\m.t->T u t = in L 2 , it is a potential and the 
trajectories have null limits at T. Therefore, this approximation also holds 
uniformly on the trajectories, on the closed interval [0, T], 

lim sup \u?{W t )-u t (W)\=0 P m -a.s. 

n ^°°0<t<T 

The function u n solves the equation (dt + L)u n + f n = with the condition 
11^ = and its backward representation is 

rT d rT 



v?(W t )= [ fs(W s )ds-J2 [ d lU n s (W s )dWl. 
Jt i=1 Jt 

If we set A™ = J * f™(W s )ds, after conditioning, this representation gives 
(*) unW t )=A%-A?-J2 f ' d iU ^W s )dWi=M m [A^/T t ]-A^ 

■ -i J t 



t=i 

In particular, one deduces 

d r x 



u n (W Q ) = E m [A%/F ] =AV t -Y j f d iU n s (Xs) dWl 

i=i Jo 

Also from the relation (*), it follows that 

E m (yl™ - A™) 2 = E m I u?(W t ) + > ' / diU n s {W s ) dWl 




T \ 2 



(**) 



\uf\\^ + / WVvZWUs. 



H 112 

A similar relation holds for differences, in particular one has 

E m {A n _ ^2 = ||u n _ u kj2 + 2 /" T y y(< _ u *)||2 ds> 

J 
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On the other hand, the preceding lemma ensures that hin^^oo aU a = I, 
in the space L 2 ([0,T[;H 1 (R d )), which implies 



lim f \\V(u? -u t )\\ 2 2 dt = 0. 



n^O 







These last relations imply that there exists a limit lim n A T =: At in the 
sense of L 2 (F m ). 

Let us denote by M n = (M t n ) tg [ 0i r],M = (M t ) t £[o,T} the martingales given 
by the conditional expectations M? = W n [A T /T t ]',M t = W l [A T /T t ]. Then 
one has lim,,,-^ M n = M, in L 2 (P m ), and hence 

lim E m sup |M t n - M t \ 2 = 0. 

n->oo Q< t <T 

Then the relation uf(Wt) = M™ — Af shows that the processes A n ,n £ N*, 
also converge uniformly on the trajectories to a continuous process A = 
(A t )te[o,T]- The inequality 

sup \A2-A t \<A T + \A T -A T \ 

0<t<T 

ensures the conditions to pass to the limit and get 

lim E m sup \A?-A t \ 2 = 0. 

n-Kx> Q<t<T 

Passing to the limit in the relations (*) and (**) one deduces the relations 
(i), (ii) and (hi). 

In order to check the relation (iv) from the statement, we observe that 
the relation is fulfilled by the functions u n , 

cT / -i \ rT 



Lp(s, x)f n (s, x) ds dx 

T 



= E m f V {s,W s )dA n s 
Jo 

where (p is arbitrary in T>t- In order to get the relation (iv), it would suffice 
to pass to the limit with n — > oo in this relation. The only term which poses 
problems is the last one. The uniform convergence on the trajectories implies 
that, P m -a.s., the measures dAf weakly converge to dA%. Therefore, one has 

lim / <p t (W t )dA?= [ n(Wt)dA t P m -a.s. 
n ^°°Jo Jo 



On the other hand, one has 
rT 

.2m^ , a2 



MW t )dA? 



o 



< sup rt{W t ) + A T + \A T -A T \\ 

0<t<T 
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By Ito's formula and Doob's inequality, one has 

/ r-T \ 2 

E m ( sup ^ 2 (t,W t )) <4\\<po\\ 2 + W m ( / \dt<p(t,W t )\dt 

MKt<T y \Jo 

+ 16E m f \Vip\ 2 (t,W t )dt 
Jo 

+ 2E m (£ \Acp\(t,W t )dt 



<4\\<p \\ 2 + 4T [ \\d m \\ 2 2 dt + 16 [ \\V<p t §dt 
Jo Jo 



+ 2T [ \\A(pt\\l<tt<oo. 
Jo 



The preceding estimate ensures the possibility of passing to the limit and 
deducing that 

limE m f ip(s,W s )dA^ = E m [ <p(s,W s )dA s , 
n Jo Jo 

and thus we obtain the relation (iv). 

Let us now consider the converse. Assume that uG F and A is a continu- 
ous increasing process adapted to (J~t)t€[o,T] an d satisfying the relation (i). 
In order to simplify the subsequent notation, it is convenient to extend our 
given function by putting m = for t>T. Now, we shall show that 

(7) P r (u t+r )<u t , te[0,T],r>0. 
By the Markov property, one gets 

P r u t+r (W t ) = E w *[ Ut+r (Wr)) = E m [u t+r (W r+t )\T t ] 

= E m [E m [A T \F t+r \ - A t+r \F t \=E m [A T \F t ] - A t+r , 

where the last line comes from the relation (i). This shows that 

P r u t+r (W t )< Ut (W t ) P m -a.s. 

and as the distribution of Wt under P m is m, we deduce the inequality (7). 
Moreover, this inequality shows by iteration that if r < r', then 

(8) Pr'Uf+r' < P r U t +r- 

By the properties of the semigroup density and since t—>ut is continuous 
with values in L 2 , it follows that, for each r > 0, P r ut+ r , t E [0,T], has a 
continuous version in [0, T] x R d defined by 

u r (t,x)= q r (x,y)u t+r (y)dy. 
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The inequality (8) shows in fact that is supermedian with respect to 
(Pt)t>o and, because of continuity, in fact it is excessive. Then u = lim r _>.o^ r 
is also excessive and since lim rH >o P r ut+ r = 1H, in L 2 , clearly u is a version of 
u. The process (ut(Wt))t£\o,T] 1S a ccu g supermartingale, and more precisely 
a potential. By the relation (i), this process admits a continuous version. It 
follows that itself is continuous and, as a consequence, one has the following 
convergence, uniformly on the trajectories: 

lim sup \u r t (Wt)-u t (Wt)\ = F m -a.s. 

r ^°0<t<T 

On the other hand, by the representation (i) one has 

E m sup \u t {W t )\ 2 <oo, 

0<t<T 

which leads to 

lirnE™ sup \lf t {W t )-ut(W)\ 2 = 0. 

0<t<T 

This relation implies that u is quasicontinuous, and hence it is a regular 
potential, completing the proof. □ 

It is known in the probabilistic potential theory that the regular potentials 
are associated to continous additive functionals (see [4], Section IV. 3 or [10], 
Theorem 5.4.2). In the above theorem, the additive aspect is not evident. 
In fact, it is hidden in the relation (i) of Theorem 2. This relation implies 
that, for t < s, A s — At is measurable with respect to the completion of 
o~(W r /r £ [t, s]). This can directly be proven but it also follows from the 
approximation of A by A n . For the processes A n ,n £ N, this measurability 
property obviously holds. And this measurability ensures the fact that A 
corresponds to an additive functional for the time-space process, which we 
are not explicitly using. 

The measure v from the theorem, expressed in the relation (v), is also 
completely determined by the relation (iv), because the test functions are 
dense in C c ([0,T] x IR d ). A natural question now is whether one Radon mea- 
sure on [0, T] x M. d can be associated via the relation (iv) from the theorem 
to two distinct potentials. The answer is that there is only one such potential 
and more precisely it can be directly expressed with the density qt(x,y) in 
terms of the measure, as one can see from the next lemma. 

Lemma 3. Let u be a regular potential and v a Radon measure on [0, T] x 
W 1 such that relation (iv) holds. Then one has 

(cj),ut)= / / ( / <j){x)q s - t {x -y)dx)v{dsdy) 



for each (j) G L 2 (M d ) and t G [0,T]. 
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Proof. We first remark that the relation (iv) is in fact equivalent to 
the following more explicit one 

(ut,(p t ) + J (^(Vu s ,V<p s ) + (u s ,d s (p s )\ ds = J t J a ^(s,x)u(dsdx), 

with any tp £ T>t and t 6 [0, T\. 

Clearly, it is sufficient to prove the lemma for (f> G C c (R d ) such that (p>0. 
Then we set ip(s, y) = J Rd cj)(x)q s -t(x — y) dx, for s € [t, T] and y € M. d . Then 
V> s = P s -t4> and the map s -> ip s is in C 1 (]t,T];L 2 (R d )) and d s ^ = \M>s- 
Let ry € C C (1R+) be a decreasing function such that rj = 1 on the interval [0, 1] 
and 7] = for x > 2. Set r/ n (x) = ??(^ )> so that (f? n )neN is an increasing 
sequence in C c (M. d ) with limit l K< j. For each fixed n, the function r/ n ^ can 
be approximated by convolution with smooth functions and then by test 
functions from T>t, and consequently we may write the relation (iv) in the 
form 

(ut,T) n ipt) + Q(Vu s , V(rj n tp s )) + (u s ,r) n d s tp s )j ds 

T 

7] n (x)ip(s, x)v(ds dx). 

Then it is easy to see that we may pass to the limit with n — > oo, in this 
relation too. Then we get 

(u t ,ip t ) + J^ Q(Vn s ,VV's) + (u s ,d s i/j s )\ ds = J t J ^(s,x)u(dsdx), 
which becomes the relation asserted by the lemma, on account of the relation 

We now introduce the class of measures which intervene in the notion of 
solution to the obstacle problem. 

Definition 2. A nonnegative Radon measure v defined in [0,T] x M. d 
is called regular provided that there exists a regular potential u such that 
the relation (iv) from the above theorem is satisfied. 

As a consequence of the preceding lemma, we see that the regular mea- 
sures are always represented as in the relation (v) of the theorem, with a 
certain increasing process. We also note the following properties of a regular 
measure, with the notation from the theorem. 

1. A set B £ B{[0,T] x R d ) satisfies the relation u(B) = if and only if 
ffl B {t,W t )dA t = 0¥ m - a .s. 
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2. If a set B G B(]0,T[xR d ) is polar, in the sense that 

¥ m ({u}£n'\3t£[0,T],(t,W t (oj))£B}) = 0, 

then v(B) = 0. 

Borcl measurable and such that > 
i/) 2 (t, x), dt <g> dx-a.e., and the processes (ipt(^t))t£[o,T] >* = 1; 2, are a.s. 
continuous, then one has ^(t/; 1 < ■0 2 ) = 0. 

2.3. Hypotheses. Let i? = (B t )t>o be a standard (^-dimensional Brown- 
ian motion on a probability space {Q,T B ,F). So = (Bj,...,B? ) takes 
values in R rf . Over the time interval [0, T] we define the backward filtration 
(J r ^ r ) sg [ 0j T] where Ff T is the completion in F B of <r(i? r — l? s ; s <r<T). 

We denote by %t the space of i/ 1 (R rf )- valued predictable and Tf T - 

adapted processes (ut)o<t<T such that the trajectories t—¥ Ut are in F a.s. 
and 

\\u\\\ < CO. 

In the remainder of this paper, we assume that the final condition is a 
given function in L 2 (R rf ) and the functions appearing in equation (1), 

/ : [0, T] x n x R d x R x R d -> R, 

5 = (g u ...,g d ): [0,T] xSlxK^xtxR^ R d , 

h = (/ti,...,/t d i):[0,T] x x R d x R x R d ^R d \ 

are random functions predictable with respect to the backward filtration 
(Ffir)telo,Ti- We set 

/(-,-, -,0,0):=/°, 9 ( v ,-,0,0):= 5 ° = (<?°,...,<?°) 

and 

/i(-,-,-,0,0):=/i° = (/i?,...,/ i ° 1 ) 
and assume the following hypotheses. 

Assumption (H) . There exist nonnegative constants C,a,/3 such that 

(i) \f(t,u,x,y,z) - f(t,u,x,y',z')\ < C(\y - y'\ + \z- z'\). 

(ii) (E J til^(*^^,y^)-^(*^,^2/'^ / )l 2 ) 1/2 <C'|y-y / |+/3|z-z'|. 

(iii) (Eti \9i(t, x, y, z) - gi {t,oj, x, y f , z')| 2 ) 1/2 <C\y - y'\ + a\z - z'\. 

(iv) The contraction property (as in [7] ) : a + < | . 

Assumption (HD2). 

E(||/ ||| 2 + lb |lL + ll^°llL)<oo. 
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Assumption (HO). The obstacle v(t,u,x) is a predictable random 
function with respect to the backward filtration (F^ T ). We also assume 
that 1 1—} v(t, uj, Wt) is P®P m -a.s. continuous on [0,T] and satisfies 

v(T, ■)<$(■)• 

We recall that a usual solution (nonreflected one) of the equation (1) with 
final condition ut = is a processus u G %t such that for each test function 
</3 G T>t and any Vi G [0, T] , we have a.s. 



(9) 



(u s ,d s ip s ) + ^(Vu s ,V(p s ) + (g s ,X7ip s ) 

= / (f s ,ip s )ds+ / (h s ,ip s ) ■ dB s . 
Jt Jt 



ds - ($,(pr) + (ut,y>t) 



By Theorem 8 in [7], we have existence and uniqueness of the solution. More- 
over, the solution belongs to Ut- We denote by U($>, f,g,h) this solution. 

Remark 1. Let L = . did 13 dj be an elliptic operator in divergence 
form, with the matrix a = {a lJ ) : M. d — > W 1 x M. d being symmetric, measurable 
and such that 

Aiei 2 <^^'(x)re j <A^i 2 

ij 

for any x,£ G M d . If instead of the operator in our equation (1), we 
had the operator L, then the contraction condition (iv) of Assumption (H) 

would be replaced by a + y < A (this ensures the contraction condition 
as formulated in [7]). Then the time change t — > ^t' yields a one to one 
correspondence between the solutions u of the equation 

du t + [Lut + f t {u t ,Vu t ) + div g t (u t ,Vu t )} dt + h t (u t ,Vu t ) ■ dB t = 0, 

over [0,T] and the solutions ut = «i/(2A)t satisfying the equation 

du t + [\&u t + ft(ut,Vu t ) + dwgt(ut,Vu t )]dt + ht(ut,Vu t ) ■ dB t = 0, 
over the interval [0,2 AT], with the transformed coefficients 

f(t,x,y,z) = ^f(^t,x,y,z^,h(t,x,y,z) = ( 2 \)V* h {^h t,x,y,z 
9i(t, x, y> z ) = J^ (si x, y, z) + ^2 a%3 ( x ) z j ~ Az i) , i = l,...,d, 
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and the transformed Brownian motion Bt = (2A) 1 / 2 i? 1 /( 2 A)i^ £ [0,2 AT]. 
This can be checked just by direct calculations using the above definition of a 
solution. Moreover, if one writes L in the form Lu = AAu — div^Vw), where 
7 = (t*"') i s a matrix with the entries 7 lJ (x) = A<5 U — a l i(x), i,j = 1,. . . ,d, 
then one has 

0<7 = AJ-o<(A-A)I, 

in the sense of the order induced by the cone of nonnegative definite matrices. 
This implies that one has 

|7(z)£|<(A-A)|£| 

for any x, £ G M d . Then it easy to deduce that <7t(a;, y, z) = ^(ffi/^A)* ( x > Vi z ) + 
7(x)z) fulfils condition (iii) of Assumption (H) with a constant a = -^(a + 
(A — A)). On the other hand, one can see that h satisfies condition (ii) with 
P = ( 2A )i/2 Pi so t^ a t the condition a + ^- < A, ensures a + ^- < |, which 
is condition (iv) of our Assumption (H). Therefore, we conclude that our 
framework covers the case of an equation that involves an elliptic operator 
like L, because the properties of the solution u are immediately obtained 
from those of the solution u. 

2.4. Quasi- continuity properties. In this section, we are going to prove 
the quasi-continuity of the solution of the linear equation, that is, when 
f,g,h do not depend of u and Vn. To this end, we first extend the double 
stochastic Ito's formula to our framework. We start by recalling the following 
result from [7] (stated for linear SPDE). 

Theorem 3. Let u £ %t be a solution of the equation 
du t + \ Aut dt + (ft + divg t ) dt + h t d~B t = 0, 
where f,g,h are predictable processes such that 

E / [WftWl + WdtWl + IIMl] dt < 00 and \\®\\2 <o °- 

Jo 

Then, for any < s < t < T, one has the following stochastic representation, 
F®F m -a.s., 

t rt 



(10) 



u(t, W t ) - u(s, W s ) = Y J f d lU (r, W r ) dW l r - [ f r (W r ) dr 
- i y g*dW- J h r {W r ) ■ tB r . 
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We remark that Tt and Tq T are independent under P ® P m and therefore 

in the above formula the stochastic integrals with respect to dWt and clWt 
act independently of J-q T and similarly the integral with respect to dBt acts 
independently of Ft- 

In particular, the process (ut(Wt))t£\o : T] admits a continuous version which 
we usually denote by Y = (^t)te[o,T] an d we introduce the notation Z t = 
Vut(Wt). As a consequence of this theorem, we have the following result. 

Corollary 1. Under the hypothesis of the preceding theorem, one has 
the following stochastic representation for u 2 , P<g>P m -a.e., for any <t <T, 



u 2 t {W t )-<$> 2 {W T ) = 2 J 



(11) 



u s f s (W s )-^\Vu s \ 2 (W s ) 



-(Vu s ,g s )(W s ) + -\h s \ 2 (W s ) 



ds 



+ J ( Ur 



■g r ){W r ) * dW r - 2^ ! (u r diU r )(W r )dW* 
i 

+ 2^ (u r h r )(W r )-d~B r . 
Moreover, one has the estimate 

EE m ( sup |y s | 2 ) +E / ||Vu s ||lds 

H<s<T ' [Jt 

l + ^j t [\\fs\\l + \\9s\\l + \\h s \\l]ds 



(12) 



for each t G [0, T] . 

Remark 2. With the notation introduced above, one can write the re- 
lation (11) as 

\Y t \ 2 + f \Z r \ 2 dr=\Y T \ 2 + 2 [ Y r f r (W r )dr-2 [ (Z r ,g r (W r )) dr 
Jt Jt Jt 

(13) Y r g r (W r ) * dW r - 2J2 jT Y r Z^ r dW l r 

+ 2 f Y r h r (W r )-(W r + [ \h r \ 2 (W r )dr. 
Jt Jt 
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Proof of Corollary 1. Assume first that g is uniformly bounded and 
belongs to {HtY, so that E || div g t \\ 2 dt < oo. Then we may represent the 
solution in the form 

u t (W t )-u s (W s ) = J2 f d iUr {W r )dW;- [\f r (W r ) + dwg r (W r )]dr 

J s J s 



h r {W r )-dB r . 
By Lemma 1.3 of [16], we may write 

u 2 t {W t )-u 2 s {W s ) = -2 f [u r (f r + dw g r )(W r ) - \Vu r \ 2 (W r ) - \h r \ 2 (W r )]dr 

J s 

+ 2^2 [ {u r diU r ){W r )dW l r -2 f (u r h r )(W r )-a 7 B r . 

j J s J s 

On the other hand, by Lemma 3.1 of [17], one has 

di\{u r g r )(W r ) dr = I u r g r (W r ) * dW r , 



so that the preceding relation immediately leads to the relation (11). Then 
the standard calculations of BDSDE involving Young's inequality, BDG in- 
equality and Gronwall's lemma give the estimate (12). 

Finally, to obtain the result with general g one proceeds by approximation. 

□ 



In the deterministic case, it was proven in [17] that the solution of a 
quasilinear equation has a quasicontinuous version. Here, we shall prove 
the same property for the solution of an SPDE as is stated in the next 
proposition. 



Proposition 1. Under the hypothesis of Theorem 3, there exists a func- 
tion u:[0,T] x x R rf R which is a quasicontinuous version of u, in 
the sense that for each e > 0, there exits a predictable random set D e C 
[0,T] x x M. d such that V-a.s. the section is open and u(-,uj,-) is con- 
tinuous on its complement (D^) c and 

¥®F m ((uj,uj')\3te [0,T] s.t. (t,u,Wt(uj')) G D e ) < e. 

In particular, the process (ut(Wt))t^\o,T] has continuous trajectories, P®P m - 
a.s. 
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Proof. Let us choose k € N with k > i , so that the Sobolev space 
H k (R d ) is continuously imbedded in the space of Holder continuous func- 
tions C 7 ^), with 7 = 1 + [f ] - f . We first assume that </> £ H k (R d ) and /, 
gi,...,g d ,hi,...,h d i belong to L 2 ([0,T] x Q; H k (R d )). By Theorem 8 in [7], 
applied with respect to the Hilbert space H k (R d ), one deduces that the solu- 
tion u = U(<&, f,g,h) has the trajectories t — > Ut(uj, •) continuous in H k (R d ) 
which implies that they are in C[[0,T] x R d ). On the other hand, we have 
from (12) the following general estimate 

EE m ( sup u(t,W t ) 2 ) <ce[||*|||+ [ T (\\ft\\l + \\gt\\l + \\ht\\l)dt . 

v 0<t<T ' Y JO 

Now, for general /, g, h), one chooses an approximating sequence of data 
f n , g n , h n ) which are H k (R d ) -valued and such that 

E^* n - $ n+1 n 2 + J\\\fr - fr +l \\i + m - 9t +1 \\i + m- h^'wi] 
i 

< — . 

- 2 n 

Let u n be the sequence of P-a.s. continuous solutions of the equation associ- 
ated to (<F\ f n ,g n ,h n ). Then set E e n = {\u n - u n+l \ > e} and D% = U n > fc E e n . 
Then we have 



e 2 F®F m {(uj,uj')\3t£ [0,T] s.t. (t, u, W t {J)) £ E e n ) 
< EE r 



sup {v%(W t )-u n t + \W t )) 2 

0<t<T 



C 

< . 

- 2" 



Further, one takes e = -y to get 



oo 4 

P ® P m (( W , o/)|3t G [0, T] s.t. (i, u,, W t {J)) eDl)<J2^- 

n=k 

This shows the statement. □ 



We also need the quasicontinuity of the solution associated to a random 
regular measure, as stated in the next proposition. We first give the formal 
definition of this object. 

Definition 3. We say that u £ T~Lt is a random regular potential pro- 
vided that u(-,uj, ■) has a version which is a regular potential, P(dcj)-a.s. 
The random variable v : Q — > M([0, T] x R d ) with values in the set of regular 
measures on [0,T] x R d is called a regular random measure, provided that 
there exits a random regular potential u such that the measure v(uS){dtdx) 
is associated to the regular potential u(-,oj,-), P(dw)-a.s. 
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The relation between a random measure and its associated random regular 
potential is described by the following proposition. 

Proposition 2. Let u be a random regular potential and v be the asso- 
ciated random regular measure. Let u be the excessive version of u, that is, 
u(-,uv) is a.s. an (Pt)t>o-excessive function which coincides with u(-,u,-), 
dtdx-a.e. Then we have the following properties: 

(i) For each e > 0, there exists a (J- ^ T ) t ^ ^-predictable random set 
D e C [0, T] x Q, x ]R rf such that P-a.s. the section is open and u(-,uj, •) is 
continuous on its complement (D^) c and 

F®P m ((cu,u')\3te [0,T] s.t. (t,uj,W t {uj))eD £ J<e. 

In particular, the process (ut(Wt))t^\o,T] has continuous trajectories, P(g)P m - 
a.s. 

(ii) There exists a continuous increasing process A = (-At)t£[o,Ti defined 
on f2 x Q' such that A s — At is measurable with respect to the P ® P m - 
completion of Tf T V a(W r /r £ [t, s]), for any < s < t < T, and such that 
the following relations are fulfilled a.s., with any ip £T> and t € [0,T]: 

(a) (ut,<pt) + J 4 T (|(Vn s , V(p s ) + (u s ,d s (p s ))ds = f Rd ip(s,x)i'(dsdx), 

(b) u t {W t ) = E[A T \T t V J*] - A t , 

(c) u t (W t ) = At — A t — Eti !t 9 MW S ) dWl, 

(d) Kill + It \Nu s \\ 2 2 ds = E m (A T - A t ) 2 , 

(e) v(<p)=W n fZ'<p(t,Wt)dA t . 

Proof. The proof of this proposition results from the approximation 
procedure used in the proof of Theorem 2. 

(i) Let r > 0. The process If = (ut)te[o,T]i defined by lf t = P r ut+ r , has 
the property that (t,x) —>v^ is jointly continuous P-a.s. We also have 

limEE m sup \u r t (W t ) -u t (W t )\ 2 = 0, 

0<t<T 

by the arguments used at the end of the proof of Theorem 2. This one 
concludes as in the proof of the preceding proposition. 

(ii) The construction of the increasing process described in Theorem 2 
holds globally for a random regular potential producing on a.e. trajectory 

the increasing process corresponding to u(-,u), •). □ 

We remark that, taking the expectation of the relation (i.i.d.) of this 
proposition one gets 

EE m (A T )=E^\\u \\l + ||Vut||i<ftY 
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3. Existence and uniqueness of the solution of the obstacle problem. 

3.1. The weak solution. We now precise the definition of the solution 
of our obstacle problem. We recall that the data satisfy the hypotheses of 
Section 2.3. 

Definition 4. We say that a pair (u, v) is a weak solution of the obstacle 
problem for the SPDE (1) associated to ($, /, g,h,v), if: 

(i) u € Ht and u(t, x) > v(t, x), dF dt <8> dx a.e. and u(T, x) = 
dF <g) dx a.e., 

(ii) v is a random regular measure on (0,T) x M. d , 
(hi) for each <p S Dt, and i £ [0, T], 

rTr ! 1 

[u s ,d s ip s ) + -(Vn a , Vy s ) ds - ($,<pt) + («t,¥>t) 

(14) = J [(f s (u s , Vu s ),(p s ) - {g s (u s , Vu s ), V(f s )} ds 

+ / (h s (u s , Vu s ),ip s ) ■ (iB s + / / cp s (x)u(ds,dx), 

Jt Jt J~R d 

(iv) if u is a quasicontinuous version of u, then one has 

r-T 

(u s (x) — v s (x))u(ds dx) = a.s. 

i 

We note that a given solution n can be written as a sum u = U\ + v,2, 
where u\ satisfies a linear equation u\ =U(Q, f(u,S/u),g(u,S/u),h(u,Vu)) 
with f,g,h determined by u, while U2 is the random regular potential cor- 
responding to the measure v. By Propositions 1 and 2, the conditions (ii) 
and (iii) imply that the process u always admits a quasicontinuous version, 
so that the condition (iv) makes sense. We also note that if u is a quasicon- 
tinuous version of u, then the trajectories of W do not visit the set {u < v}, 
P<g>P m -a.s. 

Here is the main result of our paper. 

Theorem 4. Assume that the Assumptions (H), (HD2) and (HO) hold. 
Then there exists a unique weak solution of the obstacle problem for the 
SPDE (1) associated to (<f>,f,g,h,v). 

In order to solve the problem, we will use the backward stochastic differ- 
ential equation technics. In fact, we shall follow the main steps of the second 
proof in [9], based on the penalization procedure. 

The uniqueness assertion of Theorem 4 results from the following com- 
parison result. 
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Theorem 5. Let &,f',v' be similar to $, /, v and let (u,u) be the so- 
lution of the obstacle problem corresponding to (<E>, f,g,h,v) and (u', is') the 
solution corresponding to (<&',f',g,h,v'). Assume that the following condi- 
tions hold: 

(i) $ < <£>', dx(g>dF-a.e. 

(ii) f(u, Vit) < f'(u, Vu), dt dx ® P-o.e. 

(iii) v < v' , dt dx ® P-a.e. 

T/ien one has u < u' , dtdx ® P-a.e. 

Proof. The proof is identical to that of the similar result of El Karoui 
et al. ([9], Theorem 4.1). 

One starts with the following version of Ito's formula, written with some 
quasicontinuous versions u,v! of the solutions u,u' in the term involving the 
regular measures v, v\ 

E||(«t - u' t ) + \\ 2 2 +K^ \\V(u s - u' s ) + \\lds 

= n($-$y\\l + 2Ej\(u s -uXJs(u s yu s )-f s (u' s ,Vu' s ))ds 
+ 2E f [ {u s -u' s ) + {x){v -v')(dsdx) 

Jt JR d 

+ 2E^ (V(u s -u' s ) + ,g s (u s ,Vu s )-g s (u' s ,Vu' s ))ds 

+ EJ \\h s (u s ,Vu s ) - h s (u' s ,Vu' s )\\ 2 2 ds. 

We remark that the inclusion {u > u'} C {u > v} U {v > v'} U {v 1 > u'} and 
the fact that the set {v > v'} U {V > u'} is not visited by W, imply that 
is(u > u') = 0, a.s. Therefore, 

I [ (u s -u' s ) + (x)(u-u')(dsdx)<0 a.s. 
Jt Ju. d 

and then one concludes the proof by Gronwall's lemma. □ 

3.2. Approximation by the penalization method. For n S N, let u n be a 
solution of the following SPDE 

dv?(x) + \ Au? (x) dt + f(t,x, u?(x),Vu?(x)) dt 

(15) + - v t {x))~ dt + div(#(i, x, u?(x), Vu i n (s))) (it 

+ fc(t, x, i#(<f), V<(x)) -(W t = 



OBSTACLE PROBLEM FOR STOCHASTIC PDE'S 



23 



with final condition = <J>. 

Now set f n (t, x, y, z) = f(t, x, y, z) + n{y — vt(x))~ and i> n (dt, dx) := n(u™(x) ■ 
vt(x))~ dtdx. Clearly for each n G N, f n is Lipschitz continuous in (y,z) uni- 
formly in (t,x) with Lipschitz coefficient C + n. For each n € N, Theorem 
8 in [7] ensures the existence and uniqueness of a weak solution u n G Ht of 
the SPDE (15) associated with the data ($, f n ,g,h). We denote by Y t n = 
u n (t, W t ), Z n = Vu n (t, W t ) and S t = v(t, W t ). We shall also assume that u n 
is quasi-continuous, so that Y n is P <8> P m -a.e. continuous. Then (Y n ,Z n ) 
solves the BSDE associated to the data (3>, f n ,g,h) 

Y t n = $(W T ) + [ T f r (W r , Y r n , Z n r ) dr + n f (Y r n - S?)~ dr 
Jt Jt 

+ \f g r (W r ,Y r n ,Z?)*dW 
+ J h r (W r ,Y r n ,Z?)-d~B r 

We define K™ = n J^iYJ 1 — S s )~ ds and establish the following lemmas. 



(16) 



Lemma 4. The triple (Y n ,Z n ,K n ) satisfies the following estimates 



EE" 1 \Y n \ 2 + A, EE'" / \Z?\ 2 dr 



< cEE m 

(17) 



\HWt)\ 2 4 / T (|/ s ° (W s )\ 2 + \g° s (Ws)\ 2 + \h° s (W s )\ 2 ) ds 



T 



+ c e EE m / \Y r n \ 2 dr + c 5 EE m ( sup | S r \ 2 

t H<r<T 
n\2 



+ 5EE m {K%- K\ 

where \ e = 1 — 2a — j3 2 — e, c e ,cs are a positive constants and e > 0, 6 > 
can be chosen small enough such that X E > 0. 

Proof. By using Ito's formula (13) for (Y n ,Z n ), we get 

\Y t n \ 2 + [ T \Z?\ 2 dr = MW T )\ 2 + 2 [ T Y s n f s (W s ,Y s n ,Z?)ds 
Jt Jt 

+ 2 f Y^dK n s -2 [ T (Z^g s (W s ,Y s n ,Z^)ds 
Jt Jt 
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(18) 



+ ^ Y?g a (W.,Yr,Z?)*dW-2Y i J t Y^Z^dW, 

+ 2j\ s n h s (W s ,Y?,ZZ)-kB s 
+ £ \h s (W s ,Y s n ,Z?\ 2 ds. 



Using Assumption (H) and taking the expectation in the above equation 
under P <g> P m , we get 



EE' m |K n | 2 +EE" 



\Z?\ 2 ds 



<E|d>(W T )| 2 + c £ EE m / [\f°(W s )\ 2 + \g° s (W s )\ 2 + \h s (W s )\ 2 }ds 



+ c £ EE m f \Y s n \ 2 ds + (2a + (3 2 + e)EE m f \Z n s \ 2 
Jt Jt 



ds 



+ -EE 7 ' 

7 



sup \S S \ 2 + 7 EE m [(K%-K?) 2 }, 

t<s<T - 1 



where e > 0, 7 > are a arbitrary constants and c E is a constant which can 
be different from line to line. We have used the inequality j7 YJ 1 dR™ > 

§7 S™ dlC™ and then we have applied Schwartz's inequality. We also have 
used the fact that under the measure P m the forward-backward integral 
J Y™g(r, W r , Y™, Z™) * dW as well the other stochastic integrals with re- 
spect to the brownian terms have null expectation under PtgiP" 1 . Finally, 
Gronwall's lemma leads to the desired inequality. □ 



Lemma 5. 



EE m [(^T - Kff) < c'[EE m \Y t n \' z + \\<S>f 2 ] 



(19) 



EE r 



\Y?\ 2 + \Z 1 }\ 2 ]ds 



+ ^l\\\fX + \\9° s ^ + \\h s \\ 2 2 ]d S 



Proof. Let now (u n ) nS N be the weak solutions of the following linear 
type equations 

du n t + I + div g t (uf, V<) dt + ht(%%, V<) • tS t = 0, 
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with final condition u% = 0. Set Y t n = u n {t, W t ) and Z n = Vu n {t, W t ). Then 
by the estimate (12), one has 



(20) 



EE r 



\y t r+ \z»\ds 



< cA, 



where A = EE m f^[\g s (W s , Y s n , Z™)\ 2 + \h s {W s , Y s n , Z™)| 2 ] ds. Since 
verifies the equation 

d t « - u?) + lA(u n - <) + / t (u?, V<) + n(u? - ut)" = 0, 

we have the stochastic representation 

yn _ yn = + ^ j^^yn ^ dr + R n _ R n 



from which one easily obtains the estimate 
EE m [(AT^ — A"™) 2 ] 

|^| 2 + |? t n | 2 + |$(W T )| 2 



< cEE 



+ f {\f s {W s )\ 2 + m 2 + \Z^)ds+ C \Z n s \ 2 ds 
Jt Jt 



Hence, using (20), we get 
EE m [(AT£ — A™) 2 ] 

<c / EE m [|y t n | 2 + |$(W T )| 2 ] 

r r T 

+ 4EE m J (\Y s n \ 2 + \Z2\ 2 )ds 



+ J [\f°(W s )\ 2 + \g° s (W s )\ 2 + \h s (W s )\ 2 ]ds 
which gives our assertion. □ 

Lemma 6. The triple (Y n ,Z n ,K n ) satisfies the following estimate 
EE m ( sup |Y s n | 2 )+EE m [ T \Z2\ 2 ds + EE m (K%) 2 

^0<s<T ' Jo 

||$|| 2 +EE m ( sup \S S \ 2 )+E [ T [\\fX + \\9° s \\l + \\h%l}ds 

HKs<T J JO 



< C 



where c > is a constant. 
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PROOF. From (17) and (19), we get 

pT 

7n\2 dr 



(1 - 5c')EE m \Y s n \ 2 + (l-2a-P 2 -e- 8d £ )EE m f \Z? 

J s 

< (1 + c'6)\\m + (c £ + 5c' £ )A + (c £ + 5c' £ )EE m f \Y r n \ 2 ds 

J s 

+ c & EE m ( sup \S r \ 2 

t<r<T 

where A = EE m jf[\f s {W s )\ 2 + \g° s (W s )\ 2 + \h° s (W s )\ 2 } ds. It then follows 
from Gronwall's lemma that 

sup EE m (\Y s n \ 2 ) + EE m ^ \Z?\ 2 dr + EE m {K^) 2 

0<s<T Js 

< cl \\\n 2 +EE m ( sup is r | 2 )+E / T [ii/ r iii+ik iil+ii^°iil]^ . 

v 0<r<T J Js 

Coming back to the equation (16) and using Bukholder-Davis-Gundy in- 
equality and the last estimates, we get our statement. □ 

In order to prove the strong convergence of the sequence (Y n ,Z n ,K n ), 
we shall need the following result. 

Lemma 7 (The essential step). 
(21) lim EE m [ sup ((Y t n - S t )~) 2 ] = 0. 

n->oo L < t < T J 

Proof. Let (ii n ) ng N be the sequence of solutions of the penalized SPDE 
defined in (15). From Lemma 6, it follows that the sequence (f(u n , Vu n ),g(u n , 
Vu n ),h(u n ,Vu n )) n€N is bounded in L 2 ([0,T] x O x E d ;M 1+d+dl ). We may 
choose then a subsequence which is weakly convergent to a system of pre- 
dictable processes (/, g, h) and, on account of the Lemma 13 in the Appendix, 
we obtain a sequence of families of coefficients of convex combinations, 
(a fc )fc 6 N, such that the sequences 

and 



h h = ^a 1 lh(u\Vu i ) 

i£l k 
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converge strongly, that is, 

limE [ T \\f t k -ft\\ 2 2 dt = 
Jo 

and similarly for g k , g and h k ,h. 

Now for i>n, we denote by u l,n the solution of the equation 

d4' n + [±Aul' n - nu l i n + ra t + f t (u\ Vm ! ) + div Vu')] (it 

(22) , 

+ /it(u\Vu l ) -dB t = 

with final condition u^ 1 = u^. By comparison (Theorem 5), we have that 
u i,n <; u i F ur ther, we set u fc = X^jg/ a k u l,nk , where n k = inf and we de- 
duce that 

(23) u k < V oju* < lim u n , 

where the last inequality comes from the monotonicity of the sequence u n . 
Moreover, we observe that u k is a solution of the equation 

(24) du k + [| Au£ - n k u k + n k v t + ft + div g k ] dt + h k ■ cLB t = 

with final condition iij, = vt- 

Now we are going to take the advantage of the fact that the equations 
satisfied by the sequence of solutions u k have strongly convergent coeffi- 
cients. Let us denote by Y k the continuous version on [0,T] of the process 
(u k (Wt))te\o } T] j f° r an Y k £ N. We will prove now that there exists a subse- 
quence such that 

(25) lim sup \Y t k -S t \=0 P®P m -a.s. 

fc^oo <t<T 

Since the equation (24) is linear, the solution decomposes as a sum of four 
terms each corresponding to one of the coefficients f k ,g k ,h k ,v. So it is 
enough to treat separately each term. 

(a) In the case where / = 0, g = 0, h = one obtains the term corre- 
sponding to v. Then the relation (25) is a direct consequence of the Lemma 
11. 

(b) In the case where v = 0, g = 0, h = 0, the representation of Y k is given 

by 

pT d, „t 

Y k = e^ nk{s - t) f k (W s )ds-y2 e- nk{s - l) diU k r {W s )dWi. 
Jt i=1 Jt 



Thus, we have 



r e -">('-<)f^w,)ds 



<^i, T «*" ! 
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This shows that lim^^ sup < t < r | jf e~ n ^ s -^ f k {W s ) ds\ = 0, P ® F m -a.s., 
on some subsequence. For the second term in the expression of Y k , we make 
an integration by parts formula to get 

/ e- nk{s - t) d i u k s (W s )dW l s =e- nk( - T - t) U^ k -U^ k + n k / U^e'^-^ds, 
Jt Jt 

where Up k = J Q S diU k {W r ) dW*. By the Corollary 3 of Section 4, we know that 



the martingales JT*'*, A; 6 N, converges to zero in L , and hence on a subse 

4 



quence we have lim/j—^oo 

sup <t<T Wt I = 0,P<g>P m -a.s. Then by Lemma 12, 



we see that for that subsequence 



lim sup 

fe^oo <t<T 



e 



- n ^ s -^diu k {w s )dwi 



-a.s. 



Therefore, the desired result (25) holds also in this case. This time we get 
lim^oo sup < t < T \Y k \ = 0, P ® P m -a.s. 

(c) In the case where / = 0, h = 0, v = 0, the representation of K fc is given 

by 

r-T d rT 

yk = / e -n k (s-t)~k * dW _X2 e -M s ~*) d^(W" s ) dWj 

Tt i=1 Jt 



i=l 



Now the proof is similar to that of the preceding case. We treat only the sec- 
ond term in the last expression. We set tj^ = ' 5* (W r ) dW™' 1 . Integration 
by parts formula gives 



^ e -n k (s-t) d <fji,k = fri,k _ e -n h (T-t)frijk _ nk f 1 fr,k e -n k (s-t) ^ 
Jt Jt 

On the other hand, the convergence g k — >■ g implies that the backward mar- 
tingale (u\' k ) te[0)T] converges to (J? gi, r {W r ) cM^) te[ o,T] in L 2 (P®P m ). 
The other terms in the above expression of Y k may be handled similarly 
by integration by parts and taking into account Corollary 4. Using again 
Lemma 12, as in the preceding case, we get the relation (25) in the form 
lim fe ^ oo sup 0<t<T \Y k \ = 0,P® P m -a.s. 
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(d) In the case where / = 0, g = 0, v = 0, the representation of Y k is given 

by 

d rT 

Y t k = - V / e - n *( s -*) Qjfi^W,) dWj + / e- nk( - s - tr h k ■ dB s . 
i=i ^ ■ / * 

On account of Lemma 10, the same arguments used in the previous cases 
work again. 

Now it is easy to see that the relation (25) holds for the general case. On 
the other hand, (23) and (25) clearly imply the relation 

lim sup (Y t n -S t )~=0 P ® P m -a.s. 

n ->°°0<t<T 

and then, since Y n is bounded in L 2 , one gets the relation of our statement. 
□ 



We have also the following result. 



Lemma 8. There exists a progressively measurable triple of processes 
(Y t ,Zt,Kt) t( z[o t T} such that 



EE" 



(26) 



sup \Y t n -Y t \ 2 + 

0<s<T 



\Z?-Z t \ 2 dt 



+ sup \K?-K t \ 2 



0<t<T 



as n — ?• oo. 



Moreover we have that (Yt, Zt, -?Q)te[o,T] satisfies Yt > St,Vi G [0,T] and 
J^(Y s -S s )dK s = 0, P®P m -a.e. 



Proof. From the monotonicity of the sequence (/ n )nGN and the com- 
parison Theorem 5, we get that u n (t,x) < u n+1 (t,x), <itdx<8>P-a.e., therefore 
one has Y t n < Y t n+1 , for all t G [0,T], P® P m -a.s. Thus, there exists a pre- 
dictable real valued process Y = (^)i S [o,T] such that Y™ f Y t , for all t G [0, T] 
a.s. and by Lemma 6 and Fatou's lemma, one gets 

EE m ( sup |lj| 2 ) < c. 
Moreover, from the dominated convergence theorem one has 

r T 

(27) EE m / \Y t n -Y t \ 2 dt — >0 asn^oo. 

J o 
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The relation (13) gives, for n>p, 



\y? -y?\ 2 +l \z?-z?\i 

= 2 - Yf) [f s (W s ,Y s n , Z?) - f s (W s ,Y p , Z p )] ds 

(28) +2j\Y s n -Yi)d(K?-KV) 

- 2 J T (Z- - Z p , g s (W s ,Y s n , Z») - g s (W s ,Y p , Z p )) ds 
+ J\y? ~ Y?)\g B (X 8 ,Y t n ,Z») - g s (W s ,Y p ,Z p )] * dW 
" 2 E [ T (Y s n -Y?)(Zl-Zl s )dWi 

+ ^j\Ys ~ Y?)[h s (W s ,Y s n ,Z2) - h s (W s ,Y p , Z p )\ ■ dB s 

+ ^ \h s (W s ,Y s n , Z n s ) - h s (W s ,Y p , Z p )\ 2 ds. 
By standard calculation, one deduces that 

EE m / \Z n s - Z p s \ 2 ds < cEE m / \Y s n - Y s p \ 2 
Jt Jt 

(29) + 4EE m J (Y s n -Ss)~ dRP 

+ 4EE m J\y p -S s )~ dK n s . 
Therefore from Lemma 7, (27) and (29) one gets 
EE m [ T \Y t n -Y t p \ 2 dt 

<») h T 

+ EE m J \Z?-Z p \ 2 dt — >0 asn,p^oo. 
J o 

The rest of the proof is the same as in El Karoui et al. ([9], pages 721-722), in 
particular we get that there exists a pair (Z, K) of progressively measurable 
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processes with values in R d x R such that 

EE r 



sup \Y t n -Y t \ 2 + [ \Z r t l -Z t \ 2 dt 

0<s<T Jo 



+ sup \K?-K t \ 2 

0<t<T 



as n — > oo. 



It is obvious that {Kt)t£\o,T\ is an increasing continuous process. On the 
other hand, since from Lemma 7 we have lim n _>. 00 EE m [sup 0<t<T ((Y^ 1 — 
S t )-) 2 ] = 0, then, P®P m -a.s., 

(31) Y t >s t vte[o,r], 

which yields that j Q 7 (Y s — S s ) dK s > 0. Finally, we also have f T (Y s — S s ) dK s = 
since on the other hand the sequences (l^ n ) n >o and (K n ) n >o converge uni- 
formly (at least for a subsequence), respectively, to Y and K and 



/ V" - S s ) dK™ = -n [ T ((Y s n - S s yf ds < 0. 
Jo Jo 



□ 



As a consequence of the last proof, we obtain the following generalization 
of the RBSDE introduced in [9]. 

Corollary 2. The limiting triple of processes (Yt, Zt, Kt)t^[o,T] * s a so- 
lution of the following reflected backward doubly stochastic differential equa- 
tion (in short RBDSDE): 

Y t = $(W T ) + J fr(W r ,Y r , Z r ) dr + K T - K t 

(32) + ^j\ r (W r ,Y r ,Z r )*dW 

+ J h r {W r ,Y r ,Z?)-d~B r -Y, J Z^.dW l r 



with Y t > St,\/t 6 [0,T], (Kt)te[o,T] is an increasing continuous process, Kq 
and 

(33) I (Y s -S s )dK s =0. 

Jo 





Proof of Theorem 4. Since 

f-T 



[ (\K-u p t \\l + \\Vu?-Vi4f 2 )dt = E m [ (\Y t n -Yf\ 2 + \Z?-Z?\ 2 )dt, 
Jo Jo 
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by the preceding lemma one deduces that the sequence (u n ) n ^ is a Cauchy 
sequence in L 2 (0 x [0,T];H 1 (M. d )) and hence has a limit u in this space. 
Also from the preceding lemma, it follows that dK™ weakly converges to 
dK t , F®r-a.e. This implies that 

lim f [ n(u n -v)-<p(t,x)dtdx = \imE m [ <p t (W t )dK? 
n Jo Jm.<i n Jo 

= I / ip(t,x)u(dtdx), 
Jo JR d 

where v is the regular measure defined by 

f ip(t,x)v(dtdx) = E m [ (p t (W t )dK t . 
o if Jo 

Writing the equation (15) in the weak form and passing to the limit one ob- 
tains the equation (14) with u and this v. The arguments we have explained 
after Definition 4 ensure that u admits a quasicontinuous version u. Then 
one deduces that (ut{Wt))t£\o,T] should coincide with (Yt) t £\o,T]> P<S>P m -a.e. 
Therefore, the inequality Yt > St implies u > v , dt <g> P ® dx-a.e. and the 
relation L (Yt — St) dKt = implies the relation (iv) of Definition 4. □ 

4. Some technical lemmas. 

Lemma 9. Let f G L 2 ([0,T] xR d ;l) and denote by (u n ) n&i the sequence 
of solutions of the equations 

{d t + ±A)u n -nu n + / = VneN, 

with final condition uVp = 0. Then we have 

(34) [ T \\V^\\ 2 2 dt<c\- [ T \\ft\\ 2 2 dt+ [ T e-^^Wftgdt ■ 
Jo l n Jo Jo 

Proof. It is well known that the solution (u n ) nS M is expressed in terms 
of the semigroup Pt by 

r-T 



% = J e- n ^P s - t fsds. 



A direct calculation shows that one has 

r T 

n J e -"( s -*) P s _ t u ° s ds = u°t- u n t , 

which leads to 

rT 

(35) u n t = e-< T -^u° t + n / e"^"*) (u t ° - P s - t u° s ) ds. 
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The function = e n<yT is a solution of the equation (pt + ^A)^ 71 — 

nllP + f = where ft = e^ T ~^ ft- Therefore, one has the following estimate 
for the gradient of the first term in the expression of u n 



(36) 



[ T e- n ^\\Vu° t \\ 2 2 dt<c [ T e- 2n ^\\f t \\ 2 2 dt 
Jo Jo 



(see Lemma 5 of [7] for details). In order to estimate the gradient of the 
second term of the expression of u n , we first remark that 



II, 



Pr-tfr dr, 



so that one has 

t 



nV 



< n 



< nc 



T-t 



T-t 



|VP r /, 



t+r\\2 



■n 



t+r\\2 



dr ds 



dr ds, 



where we have used the well-known inequality 



\VP r ip\\ 2 < 



m\2 



for cp e L . 



Then we estimate the time integral of the norm of the gradient, which is the 
expression we are interested in, 



nV 



Ms-t) {u 0_p s _ tU s)ds 



dt 



T 







T-t 



ne 



o 



1 



h+r 1 2 



dr ds 



o 



dt 



T t-s r-T fs' rT-sVs' 

Jo Jo Jo Jo 

X 



ns' II i || 
-7=ll/*+r|| a 



ne ne 
1 



ft+r' W2 dtdr ds drds 



< 



WftWidt 



\ c 
-^ne- ns ds) <- 
2 In 



\\ht2dt. 



This estimate together with (36) imply the statement (34). □ 



Obviously, the lemma implies that lim^—^oo 

Jo ll VM ?lli dt = 0. We need a 
strengthened version of this relation, which is presented in the next corollary 
whose proof is easy, so you omit it. 
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Corollary 3. Let f,f n G L 2 ([0,T] x IR d ;]R),n G N, be such that 
lim.n_s.oo Jq Wft" ~ /tlli dt = 0. Then the solutions (u n ) n ^ of the equations 

(d t + ±A)u™ - nu n + f n = 0, 
with final condition uV^, = 0, satisfy the relation lim nH>00 || Vti^Hl dt = 0. 

Corollary 4. Let g n ,g G L 2 ([0,T] x &e suc/i that 

llfff - dtWldt = 0. T/ien £/ie solutions (u n ) neN o/ £/ie equations 

(d t + | A)u n - ra" + div g n = 0, 
with final condition u^ = 0, satisfy the relation linin^oo J*^ || Vti™!! 2 , <ii = 0. 

Proof. We regularize g by setting g\ t = P e gi^ for i = 1, . . . ,d, e > 0, i G 
[0,T]. Then c/ 4 e G i^(R d ) and / e = div# e is in L 2 ([0,T] x R d ;R). Moreover, 
we have lim e _ > o Jq \\gf — gt ||| dt = 0. Let u e ' n be the solution of the equation 

(d t + ±A)u e ' n - nu^ n + f e = 0, 

with final condition uj? = 0. By Lemma 5 of [7], one has 

[ T ||V_? - V<' n || 2 < c / T ||# - c^Hictt 
Jo Jo 

f-T 



< C [ (|i_?-_tiil + n_!-_tiii)rft. 

Jo 



On the other hand, Lemma 9 implies, for e fixed, lmin^oo j Q \\ Vit^ll 2 ; dt = 0. 
From these facts, one easily concludes the proof. □ 



Lemma 10. Let h,h n ,n G N, be ~L 2 (M d ;M dl ) -valued predictable processes 
on [0,T] with respect to (F^ T ) t >Q and such that 

E ||/i t |||_t<oo, E / Wh^Wldt <oo 
Jo Jo 

and 

rT 



lim E f \\h? - h t f 2 dt = 0. 

rwoo J Q 

Let (u n ) n £f$ be the solutions of the equations 

du n t + [±Au[ l - nu T l] dt + h T l-o r B t = 0, 
with final condition uS, = 0, for each n G N. Then one has 



T 

lim / ||V_?||l_t = 0. 
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Proof. We regularize the process h by setting h\ t = P e hij for i = 1, . . . , 
di,e>0, t € [0,T]. Then h\ t E ff<}(R d ) and E r T ||V/t£||| < oo and 
lim^oE Jq ll^l — /if Hi eft = 0. Let u e >" be the solution of the equation 

du € t ,n + T;Au e t ,n — nu t ,n + h\ ■ clBt = 

with final condition u^ n = 0, for each n G N. The relation (hi) of Proposition 
6 in [7] written with respect to the Hilbert space H = HQ(R d ) takes the form 

2 rT 



E 



Ve,n\\2 i 
U t ' l| 2 + 



2 s 



ds + n 



In particular, one has 



T 



\VuT\\ 2 ds< 



n 



\VuT\\lds 



\Vhlgds. 



E 



\Vh e s \\lds. 



Now we write the relation (hi) of Proposition 6 in [7] for the solution u n — u e ' n 
with respect to the Hilbert space H = L 2 (M, d ), 



E 



I n e,ni|2 1 



|Vu?- VuT \\ 2 , ds + 



n 



T 



\uv. 



<' n \\lds 



:E 



\K-K\\ 2 2 ds. 







In particular, one obtains 
j-T 

E 



V«?- VuV n \\l ds <E 







1^-^111^. 



From this and the preceding inequality, one deduces 



lim sup E 



\Vu™\\lds<E 



\h s — hgW^ds 



Letting e — >• 0, one deduces the relation from the statement. □ 



Lemma 11. Let v: [0,T] xl d ->R be a function such that the process 
(vt(Wt))te[o,T] admits a version S = (St)te[o,T] with continuous trajectories 
on [0, T] and such that the random variable S* = sup 0<t<T St satisfies the 
condition E m [,S*] 2 < 00. Let u n be the solution of the equation 

{d t + ±A)u n -nu n + nv = 0, 

with the terminal condition Uj, = vt- Let Y n = (^TOte^T] be a continuous 
version of the process (ut(^t))t&[o,T] > f or each n G N. Then the following 
holds: 



lim E r 

n— >oo 



sup \Y t n -s t \ 2 

0<t<T 



0. 



3G 
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Proof. Let us set up = e nt u™ and observe that this function is a so- 



lution of the equation 



(d t + ±A)u n + v = 0, 



with vt = e vt and terminal condition Uj, = vt- Writing the representation 
of Theorem 3 with g = h = for u"(Wf), one obtains 

u?(t, W t ) = e- nT v T - V ( T diU r ;(W r ) dW; + n [ T e- nr v r {W r ) dr, 
i=i Jt Jt 
and this leads to the representation of our process Y n , given by 



yn = jgm 



Then one has 

\S t -Y t \<E r > 
Let us denote by 

V n = sup 

0<t<T 



n ( T - t] S T + n£ e-< r -^S r dr 



Ft 



S t - e 



S t - e 



-n(T-i) 



St — n 



-n(T-t) 



St — n 



' n ^S r dr 



- n ^S r dr 



Obviously, one has V n < 2S* . On the other hand, one has for any fixed 5 > 0, 



(37) 



V n < sup \S t -S s \ + 2e- n0 S*. 

\t-s\<8 



This follows from Lemma 12. Prom the inequality (37), one deduces that 
limn^ooV™ = 0,P m -a.s., and hence from the dominated convergence theo- 
rem, one gets lim™-^ E m [y n ] 2 = 0. Since 

\S t -Y t n \ <wr[v n \F t ], 

Doob's theorem implies the assertion of the lemma. □ 
Finally, we mention the following calculus lemma. 
Lemma 12. Let <p € C([0, 1];R) and 6 € (0,T) ; A > 0. Then one has 



-\t 



<p(t) dt + e~ xs ^{5) - <p(0) 



< sup \<p(t) - <p(0)\ 
0<t<6 



and 



A / e- A(s -*V(s)^ + e- A(T -*V(r)-^(t) 
< sup \<p(s)-ip(r)\+2e- xs \\ip\\ 00 

\s-r\<6,s>0 
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Proof. The first inequality follows from the relation A J e~ xt dt + e~ xs = 
1. In order to check the second relation, one dominates the expression of the 
left-hand side by 

rt+S 

A / e" A(s - i V(s) ds + e~ x5 V (t + 5)- <p(t) 



+ e~ xs 



Jt+s 



and then apply the first relation to dominate the first term. □ 

APPENDIX 

The next lemma is a classical result in convex analysis, known as Mazur's 
theorem (see [5], Remark 5, page 38). We state here the result with some 
notation that is useful for our proof. Let X be a Banach space and (a; n )neN 
a sequence of elements in X. We call finite family of coefficients of a convex 
combination a family a = {cti\i S 1} where / is a finite subset of N, a, > 
for each i € I and ^ie/ a « = 1- The convex combination that corresponds to 
such a family of coefficients is the point expressed in terms of our sequence 

b y Y,i£i a i x i- 

Lemma 13. Let (x n ) ng jsj be a weakly convergent sequence of elements in 
X with limit x. Then there exits a sequence (a k )k£^ of families of coefficients 
of convex combinations, a k = {a k \i £ 1^}, such that the corresponding convex 
combinations x k = ^iG/ fe a i x i> k converge strongly to x: lim^oo \\x k — 
x\\ =0. 
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